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coupling theory and the replica approach 
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PACS 64 . 70 . Q — Theory and modeling of the glass transition 
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Abstract. - We clarify the relation between the ergodicity breaking transition predicted by mode- 
coupling theory and the so-called dynamic transition predicted by the static replica approach. 
Following Franz and Parisi [Phys. Rev. Lett. 79, 2486 (1997)], we consider a system of particles 
in a metastable state characterized by non-trivial correlations with a quenched configuration. 
We show that the assumption that in a metastable state particle currents vanish leads to an 
expression for the replica off-diagonal direct correlation function in terms of a non-trivial part 
of the replica off-diagonal static four-point correlation function. A factorization approximation 
for this function results in an approximate closure for the replica off-diagonal direct correlation 
function. The replica off-diagonal Ornstein-Zernicke equation combined with this closure coincides 
with the equation for the non-ergodicity parameter derived using the mode-coupling theory. 



Introduction. — When trying to understand a new 
physical phenomenon one usually starts with a mean field 
theory and only afterwards one tries to account for fluc- 
tuation effects. The glass transition has been studied for 
many decades but theories which aspire to describe it at 
a mean field level, starting from a microscopic descrip- 
tion of an interacting many-particle system, are quite re- 
cent. First, there is mode-coupling theory, derived in early 
eighties by Gdtze and collaborators [lJU]- Second, there 
is the replica approach formulated in the late nighties by 
Mezard, Franz and Parisi [3H6]. While the latter theory 
makes references to the former, the precise relation be- 
tween them is unclear (with the exception of disordered 
spin models). Our goal in this Letter is to bridge the gap 
between these theories. 

The mode-coupling theory of the glass transition is, es- 
sentially, a self-consistent generalization of earlier mode- 
coupling approaches dealing with long- wavelength dynam- 
ics near critical points [7] and long-time tails in time- 
dependent correlation functions [8] to the dynamics of su- 
percooled fluids at all length and time scales. The funda- 
mental approximation of this theory consists in a factor- 
ization of a time-dependent four-point correlation func- 
tion evolving with the so-called projected dynamics [5] 
(or, in the case of an overdamped system, with the so- 
called irreducible dynamics ;9|). The somewhat mysteri- 



ous nature of the four-point function involved, the uncon- 
trolled factorization approximation, and the difficulty to 
calculate corrections have made the mode-coupling the- 
ory quite controversial. On the other hand, this theory 
makes a number of predictions \2\ using as the only input 
the static structure of the system. Thus, it is verifiable 
and, therefore, it inspired a number of experimental and 
simulational studies. The major conclusion reached in pi- 
oneering simulations of Hansen and collaborators [10] is 
that the ergodicity breaking transition predicted by the 
mode-coupling theory corresponds to a dynamic crossover 
(originally named the kinetic glass transition) occurring at 
a density smaller than (or at a temperature larger than) 
that of the laboratory glass transition. Later simulations 
tested, inter alia, the wave-vector dependence of the so- 
called non-ergodicity parameter |llj and found a remark- 
able agreement with the predictions of the mode-coupling 
theory [12] . 

In contrast to the mode-coupling theory, the replica ap- 
proach, as formulated by Mezard, Franz and Parisi [3 6 
and further developed by Parisi and Zamponi P21Q3], is 
a purely static theory. It is rooted in the analysis of dis- 
ordered spin models with long-range interactions (solv- 
able exactly within mean-field theory) which exhibit the 
so-called one-step replica symmetry breaking. As showed 
in the remarkable series of papers by Kirkpatrick, Thiru- 
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malai and Wolynes [15], these models have two transi- 
tions, a dynamic transition at which time-dependent cor- 
relation functions stop decaying and a static transition at 
which the co-called configurational entropy (the logarithm 
of the number of metastable glassy states) vanishes. The 
dynamics of these models is described exactly by equa- 
tions identical to the so-called schematic mode-coupling 
equations. The dynamic transition coincides with the er- 
godicity breaking transition predicted by the latter equa- 
tions. The replica theory of Mezard, Parisi and collab- 
orators builds upon Kirkpatrick et al.'s contribution and 
deals with the thermodynamic properties of supercooled 
many-particle systems without quenched disorder. It also 
predicts the existence of two transitions: the so-called dy- 
namic transition at which the phase space of the system 
becomes split into a collection of metastable glassy states 
and the so-called ideal glass transition at which the config- 
urational entropy vanishes. The former transition is iden- 
tified with the ergodicity-breaking transition predicted by 
the mode-coupling theory. 

It should be recalled at this point that the mode- 
coupling theory was not proposed as a mean-field theory 
of the glass transition. In fact, this theory descends from 
approaches that were derived to describe distinctly non- 
mean- field features of many-particle dynamics [HIE]- 

In addition, there is another somewhat unclear aspect 
of the identification of the ergodicity breaking transition 
of the mode-coupling theory and the dynamic transition 
of the replica approach: both the locations of these transi- 
tions and the frozen-in correlations associated with them 
are markedly different. In particular, while the non- 
ergodicity parameter predicted by the mode-coupling the- 
ory agrees well with simulations [T^] , the one inferred [T3] 
from the replica approach does not 3 ,14]. Recently, it has 
been showed that the difference between the locations of 
the two transitions increases strongly with increasing spa- 
tial dimension [T6 l fT7 ] . This fact, together with a surpris- 
ing behavior of non-ergodicity parameters in high spatial 
dimensions, prompted Ikeda and Miyazaki |17j to question 
the validity of the mode-coupling theory. 

The goal of this Letter is to clarify the relationship be- 
tween the mode-coupling theory and the replica approach. 
We start with recalling that the replica approach is a gen- 
eral scheme that, in order to make quantitative predic- 
tions, has to be accompanied by specific approximations 
of the liquid state theory. Early versions of the replica 
approach [3l[5] used the so-called replicated hyper-netted 
chain approximation. Later versions [6j[13] used a molec- 
ular version of this approximation combined with the so- 
called small cage radius expansion. Finally, the most re- 
cent version [14] used an effective potential approach com- 
bined with the same small cage radius expansion. It was 
argued [6lfl4] that the latter two approaches are well suited 
to describe dense glassy systems and the ideal glass tran- 
sition. However, they are not appropriate to describe the 
dynamic transition, except in the limit of large spatial 
dimension |14j . Here we propose an alternative to the 



hyper-netted chain approximation used in Refs. [3J[5] that 
predicts the dynamic transition that coincides with the 
mode-coupling transition |18) . 

Replica approach. — We follow Franz and Parisi 
4,5 : We consider a system of particles which is tightly 
pinned to a quenched configuration. We imagine that the 
pinning potential is gradually relaxed and that at each 
step of this process the system is allowed to equilibrate 
(according to the Boltzmann measure including the pin- 
ning potential). The main assumption is that, if the den- 
sity is high enough (or if the temperature is low enough) 
at a certain point the system gets stuck in a metastable 
state characterized by non-trivial correlations with the 
quenched configuration. We assume that this metastable 
state survives in the limit of vanishing pinning potential. 
Finally, we assume self-averaging and we average over a 
distribution of quenched configurations using the equilib- 
rium Boltzmann measure. To perform this average we 
replicate the system and at the end of the analysis we 
take the limit of the number of replicas of the system r 
going to zero. Including the quenched configuration, we 
have m — r + 1 replicas and we take the limit m — > 1 . 

Ornstein-Zernicke (OZ) equations. In a 

metastable state we can use the standard apparatus of 
equilibrium statistical mechanics, including replicated OZ 
equations, 

<5n Q/3 (ri,r 2 ) = n 2 c a p (n, r 2 ) (1) 

+n^2 / dr 3 c Q7 (ri,r 3 )(5n 7( 3(r3,r2). 
7 J 

Here a, /3 and 7 are replica indices, 5n a p{ri, r 2 ) is the 
nontrivial part of the two-particle density, 5n a p{v\, r 2 ) = 
n a j3{Ti, r 2 ) — n 2 , n is the single-particle density (assumed 
to be uniform and the same in all replicas) and c a p(ri, r 2 ) 
is the direct correlation function. Note that Eqs. ([I} 
essentially define c Q( 3's and additional assumptions are 
needed in order to make use of them. 

At this point we assume replica symmetry. We recall 
that while in Ref. [5] replica non-symmetric correlation 
functions were allowed, it was found that in the metastable 
state correlations are replica symmetric. We consider only 
the metastable state and, thus, we restrict ourselves to 
replica symmetric correlation functions. We comment fur- 
ther on this issue at the end of this Letter. Specifically, 
we assume that all replica diagonal correlations are equal, 
Stioq = Snu = ... = n 2 h and coo = en = ... = c, as are 
all replica off-diagonal ones, Suqi = 5n\2 = ... = n 2 h and 
C01 = C12 = ... = c. Using the replica symmetry we eas- 
ily get the standard OZ equation for the replica diagonal 
correlations, 

h(r 1 ,r 2 ) = c(ri,r 2 ) + n J dr 3 c(ri, r 3 )h(r 3 , r 2 ). (2) 

Physically, this means that we have equilibrium correla- 
tions within each replica. 
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To get an equation for h we first re-write the replicated 
OZ equation for a =/= f3 in the following form: 



dr 3 (S(r 13 ) - nc aa (v 1 ,v 3 ))8n aB {r 3 ,T 2 ) 



= nj dr 3 c Q ^(ri,r 3 )(?ic>(r 32 ) + Jn^(r 3 , r 2 )) 

+n ^2 J dr 3 c ai (ri,r 3 )Sn 7 p(r 3 ,r 2 ) (3) 

In the m — > 1 limit from Eq. ([3]) we obtain 

rfr 3 (5(ri 3 ) - nc(r 1 ,r 3 ))h(r 3 ,r 2 ) = c(ri, r 2 ) (4) 

i J dr 3 c(r 1 ,r 3 )h(r 3 ,r 2 ) - n J dr 3 c{r 1 , r 3 )h(r 3 , r 2 ). 

We note in passing that if we define a new replica off- 
diagonal correlation function h through the following 
equation 

fe(ri,r 2 ) = y dr 3 dr 4 ((5(r 13 ) + n/i(ri, r 3 ))fr(r 3 ,r 4 ) 

x(<5(r 42 ) + nft(r 4 ,r 2 )), (5) 

we can derive from Eq. Q a somewhat simpler equation 
for h: 

h(ri,r 2 ) = c(ri,r 2 ) (6) 
-n / dr 3 dr 4 c(ri,r 3 )(5(r 34 ) +n/i(r 3 ,r 4 ))ft(r 4 ,r 2 ). 



Eqs. (06]) imply that replica off-diagonal correlation func- 
tion h can be written as a series in which the kth term is 
a convolution of the form (— l)( fc+1 )((<5 + nh)c) k (5 + nh). 

At this point we invoke the physically motivated rela- 
tion [14] between the Fourier transform of the replica off- 
diagonal correlation function, h(k), and the non-ergodicity 
parameter [11] of the mode-coupling theory, f(k), 



nh(k) = S(k)f(k). 



(J) 



Using Eq. we can re-write the Fourier transform of 
Eq. g]) as 



f(k) 
!-/(*) 



nS(k)c(k). 



(8) 



We note the striking resemblance of Eq. ([5]) and the mode- 
coupling theory's equation for the non-ergodicity param- 
eter (see, e.g, Eq. (2.17) of Ref. pQ) [19]. Specifically, 
Eq. ([5]) coincides with the equation for the non-ergodicity 
parameter iff the right-hand-side of this equation can be 
related to the long-time limit of the appropriate memory 
function. 

We would like to emphasize that the derivation of Eq. 
((8]) did not involve any approximation. The only assump- 
tions used were the existence of a metastable state and 
the replica symmetry. We believe that Eq. ([8]) is the first 
result suggesting a connection between off-diagonal direct 
correlation function in replica approach and the long-time 
limit of memory function introduced in dynamic approach. 



Metastable state. — The basic assumption of the 
replica approach is that if the density is high enough (or 
if the temperature is low enough) the system gets stuck in 
a metastable state. In the original papers this state was 
defined as a local minimum of a free energy considered as 
a functional of the replica off-diagonal correlations [3H5] 
or a function of the cage radius [5|13U14j . Here we assume 
that the metastable state is characterized by vanishing 
particle currents. We note that this condition is qualita- 
tively different than that used in earlier works in that it 
implicitly depends on the system's dynamics. This is im- 
portant since we assume metastability only with respect to 
local dynamics and we admit a possibility that non-local 
[20] or cluster [21] moves can still equilibrate the system. 
For concreteness, we will here assume that the system is 
evolving with Brownian dynamics [22 . In this case we 
will see that the condition of vanishing particle currents is 
equivalent to the assumption that the replica off-diagonal 
correlation functions satisfy the first few equations of the 
Yvon-Born-Green (YBG) hierarchy |23j . 

To motivate the assumption of vanishing particle cur- 
rents we note that the time-independence of the joint den- 
sity in replicas a and ft a / ft results in the following 
identity (to simplify notation we assume hereafter that the 
diffusion coefficient of an isolated particle has been made 
equal to 1 and that the inverse temperature 1/ksT has 
been incorporated into the inter-particle force F): 

= 9 t n a/ 3(ri,r 2 ;t) 

= -9 r Ja,/3(ri,r 2 ;t) -9 r2 j/3,a(ri,r 2 ;t) (9) 

where j a g is the joint current density in replica a and 
particle density in replica ft 



jo,/9(ri,r 3 ) = -d ri n a0 {T X ,T 3 ) 



(10) 



dr 2 F(ri 2 )n aQ/3 (ri, r 2 , r 3 ) 



and, correspondingly, is the joint current density in 
replica (3 and particle density in replica a. In Eq. (ITU1) 
n a ap is the joint two-particle density in replica a and 
single-particle density in replica ft 

The main fundamental assumption of the present ap- 
proach is that not only dtn a p — (which is required to 
have a stationary state) but moreover j Q ( g = = ]p, a - 

The vanishing of the joint current density in replica a 
and particle density in replica f3 leads to the following 
equation relating n a p and n aa Q\ 

- d ri n aa (ri,r 3 ) + J dr 2 F(r 12 )n aa 8(r 1 , r 2 , r 3 ) = 0. (11) 

In a very similar way, starting from the vanishing of the 
joint current density in replica j3 and two-particle density 
in replica a we obtain an equation relating n aa 8 and the 
joint two-particle density in replicas a and ft n aa gQ: 

j/3,oia(r 3 ,ri,r 2 ) = -9 r3 ^aa^(ri,r 2 ,r 3 ) (12) 
+ J dr i F(r 3i )n aaP/3 (r 1 , r 2 , r 3 , r 4 ) = 0. 
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We should emphasize here that equations (|1HI12[) are 
valid only for a ^ /3 and, moreover, in the limit of vanish- 
ing pinning potential. Furthermore, we should note that 
in this limit Eqs. (|imi2p are identical to the YBG equa- 
tions for n a p and n aaj 3, a ^ j3. 

Combining equations (|1HI12[) we get 

9 ri 9 r3 n 2 /i(ri,r3) = d ri d T3 n a p(rx, r 3 ) (13) 
= J dr 2 F(r 12 ) J dr 4 F(r3 4 )ra QQ/9l g(ri, r 2 , r 3 , r 4 ). 

The physical meaning of Eq. (IT31 is that correlations be- 
tween particle densities in different replicas have to be 
accompanied by correlations of forces in these replicas. 

It would be tempting to factorize the four-point func- 
tion at the right-hand-side of Eq. (Til?]) . However, due 
to non-vanishing three-particle correlations, two-particle 
densities in replicas a and /? can be correlated via single- 
particle density correlations. To take care of this, we use 
the following decomposition of ii aa ^ which is inspired by 
similar decompositions of equilibrium densities |24j : 

see Eq. (TT51) 

The first term at the right-hand-side of Eq. (fT5")) does 
not contribute to the force correlations. The second term 
describes all correlations between two-particle densities in 
replicas a and f3 mediated by single-particle correlations. 
In this term, the brackets (...) denote averaging over the 
replicated state, (Sn aa Snj) is the correlation function of a 
fluctuation of a microscopic two-particle density in replica 
a, 

Sn aa (r 1 ,r 2 ) = ^[5(r x - r ia )5(r 2 - r ia ) 

- (5(ri - r ia )6(r 2 - r ja ))] (16) 

and a fluctuation of a microscopic single-particle density 
in replica 7, 

5n T ( ri ) = 5}5(n - r. i7 ) - («5(n - r i7 )>], (17) 

i 

and {Sn 7 5ns) 1 is the inverse of the correlation of single- 
particle density fluctuations, 

((5n 7 (r 5 )(5n 5 (r 6 )) _1 = S i sS<yT ° — — - c 7 <5(r 5 , r 6 ). (18) 

n 

Finally, the third term at the right-hand-side of Eq. (TT5"]) , 
n aa8pi contains correlations between two- particle densi- 
ties in replicas a and (3 that are mediated by at least two 
single-particle correlations. We should emphasize that Eq. 
(fl"5"]) is essentially a definition of n^ a g B . 

Substituting decomposition ([15]) into Eq. (JT3J) and then 
using Eq. (1X1]) . the first equation of the standard Yvon- 
Born-Green hierarchy (23] , and the replicated OZ equa- 
tions (p} we obtain the following relation for the replica 



off-diagonal direct correlation function: 

<9 ri <9 r3 n 2 c(ri,r 3 ) = d ri d r3 n 2 c aj3 {ri, r 3 ) (19) 
= J dr 2 F(r 12 ) J dr 4 F(r 3 4)n^ (3;3 (ri, r 2 , r 3 , r 4 ) 

The only assumption used so far is the existence of a 
metastable state characterized by vanishing currents. 

Closure. — Now we are in a position to introduce the 
main approximation: a factorization of the four-point cor- 
relation function n^ a g^. For typical interactions (which 
are strongly repulsive at short distances) a straightforward 
factorization would result in divergences. To arrive at an 
alternative factorization approximation we note that Eqs. 
([5]|6]) suggest that replica off-diagonal correlations contain, 
on each side, exact single replica correlations. Thus, to 
proceed we factor out the single-replica analogs of n l " aB p 
and factorize only its remaining (middle) part: 

see Eq. (pU) 

Here n 22 is the single replica analog of n^ a8g [25], (n 22 ) _1 
denotes the inverse of n 22 and the brackets [...]f ac t denote 
factorization of the expression in the brackets, e.g., 

[ rl aa/3/3( r l) r 2,I"3,r4)]fact = Tl a p {riz)n a p (r 24 ) + (1 2). 

Explicitly, factorization approximation in Eq. (|2ip leads 
to the following expression for n^go: 

see Eq. ([24]) 

We now substitute ([24"]) into Eq. ([T9]) . take the Fourier 
transform, and follow the procedure used in the derivation 
of the so-called vertices of the mode-coupling theory [S] to 
obtain the following (approximate) identity 

kkclAO = y/^<Kk- qi -q 2 ) (25) 

x [qic(<7i) + Qacfe)] [qic(<7i) + q2c(q 2 )] ■ 

We note that the bare forces in Eq. (fl~9| got replaced by 
renormalized forces (essentially, gradients of the replica- 
diagonal direct correlation functions). Technically, this 
results from including njfl's at both "ends" of the approx- 
imate expression ([24]) . 

Next, we take the longitudinal part of Eq. ([25]) . use the 
identification ([7]) and in this way we derive an expression 
for the replica off-diagonal direct correlation function c in 
terms of the non-ergodicity parameter /, 

*i-s/w*-«-i (26 > 

xf(qi)f(q2)S( qi )S(q 2 ) (k ■ [q lC ( qi ) + q 2 c( (?2 )]) 2 

Finally, we substitute expression ([26]) into Eq. §8§ . The re- 
sulting self-consistent equation for the non-ergodicity pa- 
rameter is identical to the equation derived using the stan- 
dard mode-coupling theory. In particular, this equation 
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dr 5 dr 6 (Sn aa (r 1 ,r 2 )Sn 7 (r 5 )) (Sn 7 (r 5 )5n s {r 6 )) (Sn s (r 6 )5npf3(r 3 , r 4 )) + n™ T a/3l3 (ri, r 2) r 3 , r 4 ) (15) 



< r a/3/3( r i J r 2,r 3 ,r 4 ) 



/ dr 5 ...dri2n^'(ri,r2,r 5 ,r 6 ) 



(21) 



[n 2 r 2 r (r 5 , r 6 , r 7 , r 8 ) 1 n£* ^(r 7 , r 8 , r 9 , rio)n 22 (r 9 , r 10 ,r u ,r 12 ) x ] f n 2 r 2 r (rn, r 12 , r 3 , r 4 



implies the following relation between c and the long-time 
limit of the irreducible memory function M lrr , 

c(k) = lim M' m {k;t). (27) 

nk t->oo 

Discussion. — We have showed that the mode- 
coupling theory equation for the non-ergodicity parameter 
can be obtained by combining replica approach, condi- 
tion of vanishing currents and the factorization approx- 
imation for the non-trivial part of the replica-off diag- 
onal static four-point correlation function. One cannot 
but conclude that, at least as far as the calculation of the 
non-ergodicity parameter is concerned, the mode-coupling 
theory is equivalent to the (admittedly mean-field) replica 
approach. One could even say that in spirit it resembles 
the Kirkwood-Monroe theory of freezing 26 . 

Two features of the present development should be em- 
phasized. First, the consistency of replica approach's 
equation for the off-diagonal correlations and the mode- 
coupling theory equation for the non-ergodicity parame- 
ter is the strongest connection one could expect between 
the replica approach and the mode-coupling theory. The 
reason is that the replica approach is an inherently static 
one and the bulk of mode-coupling theory's predictions are 
concerned with dynamics. Second, since both replica ap- 
proach and mode-coupling theory share mean-field char- 
acter, it is natural to expect that both of them will in- 
volve some factorization approximations. The advantage 
of the replica derivation presented here is that it relies 
upon a factorization approximation for static correlation 
functions, albeit for a replicated system. 

There are several aspects of our approach that could be 
investigated further. First, a more accurate factorization 
approximation should be attempted. It can be argued 
that Eq. overestimates n" r aj3/3 and this leads to an 

underestimation of the mode-coupling transition density 
(or overestimation of the mode-coupling transition tem- 
perature). Second, the limit of large spatial dimensions 
should be investigated. In this limit static correlations 
simplify and it is possible that factorization approxima- 
tions could be avoided. Third, replica symmetry should 
be inspected. In particular, one should investigate break- 
ing symmetry between the Oth replica and the remaining 



r replicas. A preliminary study suggests that, at least for 
the hard-sphere interaction and with Percus-Yevick clo- 
sure for the replica diagonal correlations, this symmetry 
is not broken. Fourth, dependence of the final equation for 
the ergodicity breaking parameter on the dynamics should 
be investigated. In particular, one could try to derive the 
equation for the non-ergodicity parameter corresponding 
to Monte Carlo dynamics and to check whether the stan- 
dard equation for the non-ergodicity parameter changes. 

We thank G. Biroli and F. Zamponi for comments on 
the manuscript and gratefully acknowledge the support of 
NSF Grant CHE 0909676. 
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